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Abstract 

We report on a calculation of the 'mixed' strong and (purely) weak corrections through 
the order ag«\y to parton-parton processes in all possible channels at hadron colliders 
entering the single jet inclusive cross section. At both Tevatron and LHC, such effects 
are always negligible (below permille level) in the total integrated cross section whilst 
they become sizable in differential rates. Specifically, if such corrections are defined with 
respect to the full leading-order result of 0(ag + as^EW + a E\v)' we ^ n( ^ that, at the 
FNAL accelerator, they can reach the —5% benchmark in the jet transverse energy (at 
the kinematical limit of the machine, rendering their detection quite difficult). At the 
CERN collider, in the same observable, they exceed the —10% level already at 1 TeV and 
can reach —40% at 4 TeV, kinematic regions where such corrections will be comfortably 
observable for standard luminosity. In addition, such corrections are somewhat sensitive 
to the factorisation/renormalisation scale choice. 
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1 Introduction 

As the centre-of-mass (CM) energy of present and future hadron colliders increases well beyond 
the Electro- Weak (EW) scale, of 0(100 GeV), such as at the Tevatron at FNAL (yfs~ pP = 1.96 
TeV) and the Large Hadron Collider (LHC) at CERN (^ pp = 14 TeV), it is clear that 
the impact of higher order EW corrections, particularly of the purely Weak (W) compo- 
nent, will become more and more important phenomenologically, with respect to similar 
Quantum Chromo-Dynamics (QCD) effects [1]. The reason is twofold. On the one hand, 
the strong coupling 'constant', as, decreases with increasing energy faster than the EW one, 
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ct w = «em/ sin 2 9w (with qem the Electro-Magnetic (EM) coupling constant and 6w the weak 
angle). On the other hand, the purely weak part of higher order EW effects produces lead- 
ing corrections of the type «w log 2 (/i 2 /M^), wherein fj, represents some typical energy scale 
affecting the hard process in a given observable, e.g., the partonic CM energy V§. For large 
enough \x values, such W effects may be competitive not only with Next-to-Next-to-Leading- 
Order (NNLO) (as % rs a|) but also with NLO QCD corrections (e.g., for \i = 0.5 TeV, 
log 2 (/i 2 /M^) ~ 10). (Clearly, in actual calculations, combinatorial effects, colour /flavour fac- 
tors as well as subleading terms will have to appropriately be accounted for.) These 'double 
logs' are of Sudakov origin and are due to a lack of cancellation between virtual and real 
^-emission in higher order contributions. This is in turn a consequence of the violation of 
the Bloch-Nordsieck theorem in non-Abelian theories [2]. (See Refs. [3, 4] for comprehensive 
reviews.) 

The problem is in principle present also in QCD. In practice, however, it has no observable 
consequences, because of the final averaging of the colour degrees of freedom of partons, forced 
by confinement into colourless hadrons. This does not occur in the EW case, where the initial 
state generally has a non-Abelian charge, such as in proton- (anti) proton collisions. Besides, 
these logarithmic corrections are finite (unlike in QCD), since M w provides a physical cut-off 
for jy-emission. Hence, for typical experimental resolutions, softly and collinearly emitted 
weak bosons need not be included in the production cross section and one can restrict oneself 
to the calculation of weak effects originating from virtual corrections only [5]. By doing so, 
similar logarithmic effects, ~ «w l°g 2 (/i 2 /M§), are also generated by virtual corrections due 
to Z-bosons. Finally, in some instances, all these purely weak contributions can be isolated 
in a gauge-invariant manner from EM effects which therefore may not be included in the 
calculation (as it is the case here). (Notice that QED corrections, just like QCD ones, are not 
subject to Sudakov enhancement.) 

Furthermore, physics Beyond the Standard Model (BSM), if due to right-handed weak 
currents [6], contact interactions or compositeness [7] or new massive gauge bosons W, Z' [8, 9] 
may well manifest itself in deviations from the SM in jet quantities (which eventually emerge 
as observables of parton-parton collisions after hadronisation: see [10]- [17] for definitions and 
reviews of standard algorithms for hadronisation and jet evolution), such as jet transverse 
energy Et or di-jet invariant masses Mjj. The lower bounds on typical W, Z' masses from 
direct searches at Tevatron are model dependent but are typically around 500 — 600 GeV 
[18], whilst the LHC can access such massive gauge bosons with masses up to 4 — 4.5 TeV 
[19]. Whilst typical compositeness limits from Tevatron are very mild, compared to the EW 
scale, the LHC can probe (with planned high luminosity) jet transverse energies of about 4 
TeV or so, where the possible existence of quark substructures could be manifest. Moreover, a 
tenfold increase in luminosity, as the one currently been discussed for the so-called Super-LHC 
(SLHC) option, would give access to jets of up to E T 5 TeV [20]. 

The plan of the paper is as follows. In the next Section we describe the subprocesses 
that we have computed in terms of topologies of Feynman diagrams. Sect. 3 illustrates the 
computational techniques adopted. Sects. 4 and 5 present our results for Tevatron and LHC, 
respectively. Sect. 6 summarises our work and draws the main conclusions. Finally, in the 
appendix, we describe in detail the helicity amplitude formalism adopted here, also giving 
some specific formulae for (some of) the topologies discussed in Sect. 2. 
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2 The subprocesses 



In view of what explained in the previous section, it becomes of crucial importance to assess 
the quantitative relevance of one-loop weak effects entering via 1 0{a\a^) the fifteen possible 
2^2 partonic subprocesses responsible for jet production in hadronic collisions 2 , namely 3 : 



99 -> 99 (i) 

qq -> gg (2) 

qg -> g# (3) 

qg -> g# (4) 

99 -»■ 99 (5) 

99 -»• 99 (6) 

9Q — iQ (same generation) (7) 

qQ — > 9Q (same generation) (8) 

qQ — > gQ (different generation) (9) 

9<Q — * 9Q (different generation) (10) 

99 -»• 99 (11) 

qq — > Q<Q (same generation) (12) 

99 ~^ (different generation) (13) 

qQ — > gQ (same generation) (14) 

gQ — > gQ (different generation), (15) 



with g and Q referring to quarks of different flavours, limited to u-, d-, s-, c- and 6-type, 
all taken as massless. Whilst the first four processes (with external gluons) were already 
computed in Ref. [21], the eleven four-quark processes were tackled in Refs. [22, 23] (see also 
[24] ) 4 . Furthermore, these four-quark processes can have Infra- Red (IR) divergences, both soft 
and collinear, so that gluon bremsstrahlung effects ought to be evaluated to obtain a finite 
cross section at the considered order. In addition, for completeness, we have also included the 
non-divergent 2^3 subprocesses 



qg -> ggg (16) 

qg -> 9gg (17) 

qg — > qQQ (same generation) (18) 

qg — > qQQ (same generation). (19) 



1 Hereafter, the notation aw is used to signify purely weak effects only whilst ckew exemplifies the fact that 
both weak and EM effects are included at the given order. 

2 Note that in our treatment we identify the jets with the partons from which they originate. 

3 To clarify our notation: e.g., in the case of process (7), ll qQ — > qQ (same generation)" refers to, e.g., 
ud — > ud whereas a qQ — > qQ (different generation)" to, e.g., ub — > ub. 

4 Note that gg — ► gg does not appear through O(agaw) nor do qq' — > QQ', <j<j' — » QQ' and q<f — > QQ', if 
q' ^ q and Q'^Q. 
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(See Refs. [25, 26] for tree-level cts^EW and cc| w interference effects.) As intimated in the 
Introduction, we will instead ignore altogether the contributions of tree-level a|a:w terms 
involving the radiation of W- and Z-bosons. This presumes that the adopted jet definition 
efficiently vetoes against those gauge bosons decaying inside the jets and the detector coverage 
minimises the loss of all their decays. We will return to this matter later on in the paper. 

Before proceeding to describe the calculation we have performed, we list the 'topologies' 
involved (i.e., Feynman graphs without reference to specific couplings or internal masses) and 
briefly discuss their most salient phenomenological aspects. 

• Processes 1 (gg — > qq) and 2 (qq — > gg) (Figs. 1,2): 

The two processes with a gluon pair in the initial or final state do not exist through 
0(asaw) (as we need at least four quark-gluon vertices) thus the one-loop weak correc- 
tions to a| correspond to the leading order weak calculation. The ctgctw calculation is 
IR finite and therefore we do not include gluon bremsstrahlung diagrams. 

• Processes 3 (qg — > qg) and 4 (qg — > qg) (Figs. 3, 4): 

Comments made for the previous two processes also hold in this case. 

• Processes 5 {qq — > qq) and 6 (qq — > qq) (Fig. 5): 

The full calculation for these processes will have a contribution from order ctsaw and 
requires gluon bremsstrahlung diagrams, meaning that a subtraction procedure will have 
to be enforced in order to cancel the IR divergences between the real and virtual parts 
of the calculation (see later on for further details on the techniques adopted to achieve 
such a cancellation). Naturally, VT-exchange is here confined to self-energy and vertex 
corrections. 

• Processes 7 (qQ — > qQ (same generation)) and 8 (qQ — > qQ (same generation)) (Fig. 7): 
Hereafter, notice that we assume a flavour diagonal Cabibbo-Kobayashi-Maskawa (CKM) 
mixing matrix for quarks. The approximation is justified by the unitarity of such a ma- 
trix, the small value of the Cabibbo angle and by the fact that contributions from top 
quarks are negligible in comparison to those from light quarks. Besides, this approxi- 
mation considerably simplifies the calculation. 

• Processes 9 (qQ — > qQ (different generation)) and 10 (qQ — > qQ (different generation)) 
(Fig. 8): 

The diagrams that contribute to process 9(10) are a subset of those that contribute to 
process 7(8). As a consequences of the assumption made for the CKM matrix, we ignore 
here the interaction where a H^-emission/absorption changes a quark generation index 
and as such the W^-mediated graphs connecting the two fermion lines that contribute to 
processes 7 and 8 are ignored in 9 and 10, respectively. For the same reason, there is no 
tree-level as^w contribution here. 

• Process 11 (qq — > qq) (Figs. 9, 10, 11): 

This is possibly the most involved process of all, in terms of number of different topolo- 
gies, though W- correct ions are limited to the case of self-energies and triangles. 
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• Process 12 (qq — > (same generation)) (Fig. 12): 

The topologies are here a subset of those involved in the previous process, upon replacing 
appropriately Z- with W^-exchange between the two fermion lines. 

• Process 13 (qq — > QQ (different generation)) (Fig. 13): 

The graphs involved here are a subset of the previous process, as W^-exchange between 
the two fermion lines is not present. In addition, unlike the previous case, there is no 
allowed interference that contributes at order aga w to this process, due to the vanishing 
of the colour factors. 

• Process 14 (qQ — > qQ (same generation)) (Fig. 14): 

W-exchange between the two fermion lines here only occurs in annihilation diagrams, 
unlike all previous cases, where it takes place in scattering graphs only. 

• Process 15 (qQ — > qQ (different generation)) (Fig. 15): 

The graphs involved here are a subset of the previous process, as W -exchange between 
the two fermion lines is not present. As in process (13) there is no allowed interference 
that contributes at order asaw to this process, as all possibilities have vanishing colour 
factors. 

• Processes 16 (qg — > qqq) and 17 (qg — > qqq) (Fig. 16): 

These processes are IR finite as there are no possible bremsstrahlung diagrams we can 
write down at agCtw order. 

• Processes 18 (qg — > qQQ (same generation)) and 19 (qg — > qQQ (same generation)) 
(Fig. 16): 

The same comments as in the previous case are valid here in the context of IR finiteness. 
Note also that there are no processes analogous to processes 18 and 19 with quarks of 
different generations as there are no contributing interferences at the required order. 

3 The calculation 

Given the large number of diagrams involved in the computation, it is of paramount impor- 
tance to perform careful checks. In this respect, we should mention that our expressions have 
been calculated independently by at least two of us using FORM [27] and that some results 
have also been reproduced by another program based on FeynCalc [28]. (See the Appendix for 
more details on the procedure adopted for the tensor reduction and amplitude calculation.) 
We also find reasonable agreement with Ref. [29] in the Sudakov limit, i.e., for large invariant 
masses and transverse momenta of the final state, provided that the final state particles are 
fairly central. 

As already mentioned, IR divergences occur when the virtual or real (bremsstrahlung) 
gluon is either soft or collinear with the emitting parton and these have been dealt with by 
using the formalism of Ref. [30], whereby corresponding dipole terms are subtracted from the 
bremsstrahlung contributions in order to render the phase space integral free of IR divergences. 
The integration over the gluon phase space of these dipole terms was performed analytically 
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in d- dimensions, yielding pole terms which cancelled explicitly against the pole terms of the 
virtual graphs. There remains a divergence from the initial state collinear configuration, which 
is absorbed into the scale dependence of the PDFs and must be matched to the scale at which 
these PDFs are extracted. Through the order at which we are working, it is sufficient to take 
the LO evolution of the PDFs (and thus the one-loop running of as). 

The self-energy and vertex correction graphs contain Ultra- Violet (UV) divergences that 
have been subtracted here by using the 'modified' Dimensional Reduction (DR) scheme at the 
scale \x = M z . The use of DR, as opposed to the more usual 'modified' Minimal Subtraction 
(MS) scheme, is forced upon us by the fact that the W- and Z-bosons contain axial couplings 
which cannot be consistently treated in ordinary dimensional regularisation. Thus the values 
taken for the running as refer to the DR scheme whereas the EM coupling, oem, has been 
taken to be 1/128 at the above subtraction point. The one exception to this renormalisation 
procedure has been the case of the self-energy insertions on external fermion lines, which have 
been subtracted on mass-shell, so that the external fermion fields create or destroy particle 
states with the correct normalisation. 

The top quark entering the loops in reactions with external 6's has been assumed to have 
mass m t = 175 GeV and width r t = 1.55 GeV. The Z-mass used was M z = 91.19 GeV and was 
related to the W-mass, M w , via the SM formula M w = M z cos6 w , where sin 2 9 W = 0.232. 
(Corresponding widths were T z = 2.5 GeV and = 2.08 GeV.) Also notice that Higgs 
contributions are not included here, as required by all quarks being massless. For the strong 
coupling constant, as, we have used the one- or two- loop expression with A^- -4 -* 5 chosen to 
match the value required by the LO and NLO Parton Distribution Functions (PDFs) used. 
The latter were CTEQ6L1 (our default) plus CTEQ6L at LO and CTEQ6M at NLO [31], 
respectively. 

The fully differential cross sections for processes (1)— (19) are obtained numerically in 
FORTRAN as follows 

da = ( 27 r)V (J> - J» 1 (n (2^%) lAmxJflMixJdxidxi, (20) 

wherein = (-E 1 ^/), P?(/)) are the initial (final) particle momenta, \A\ 2 the amplitude squared 
averaged(summed) over initial (final) colours and helicities, Xi the usual Bjorken momen- 
tum fractions and /f® the proton(antiproton) PDFs. (For convenience we will refer to 
fi ( x i) f2^\ x 2)dxidx 2 - assuming a summation over all possible quark and gluon combinations 
- as the parton luminosity.) Finally, the integrations over the two- or three-body phase spaces 
and the x's have been performed using VEGAS [32]. (A simple change of kinematic variables 
yields the E T dependence of the cross section discussed below.) 



4 Tevatron phenomenology 

A detailed discussion for the case of Tevatron, geared to comparisons against existing data 
from Run 1 and 2, can be found in Ref. [23]. This is mainly a discussion on the impact of 

5 Strictly speaking, we should have amended the values taken for as in order to account for the difference 
between the MS and DR schemes, but this difference is numerically negligible. 
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(9(tts a w) one-loop virtual corrections on high Et jet samples, where an excess was initially 
found by CDF (but not DO) during Run 1 [35], with respect to the NLO QCD predictions [36]- 
[39]. It was eventually pointed out that a modification of the gluon PDFs at medium/large 
Bjorken x [40] can reconcile theory and data, even within current systematics: see, e.g., [41]. 
In fact, notice that with the most recent PDFs (e.g., CTEQ6.1M [31]), the preliminary Run 2 
data also seem to be consistent with NLO QCD, see [42] for CDF, albeit barely. The question 
as to whether 0{a\a^) effects may be required in a future comparison between theory and 
data comparisons (or in the parameterisation of the PDFs) is still open [43] and such effects 
are presently being implemented in the f astNLO package [44]. 

Here, we complement the results of [23] by studying the effects of the 0{a\a^r) one-loop 
virtual corrections over the full kinematic range of the final state jets, rather than restrict 
ourselves to the specific case of published CDF and DO data samples. Fig. 17 shows the size 
of the aforementioned corrections relative to the full LO results (as defined in the caption) for 
three sets of CTEQ PDFs [34] all taken at the factorisation/renormalisation scale // = fj,p = 
Hn = Et/2, i.e., half the jet transverse energy, the standard choice in NLO QCD simulations 
(see Refs. [36]- [39]). Weak effects are basically independent of the PDFs used and results 
are negligible for small jet transverse energy, where the differential cross section is largest 
and therefore lead to negligible corrections in the integrated one. However, the corrections in 
the differential cross sections become somewhat relevant at large jet transverse energy. For 
example, at the kinematic end Et = 800 — 900 GeV they can reach the —5% level. It is 
however doubtful whether Tevatron can reach the luminosity necessary to isolate jets at such 
large jet transverse energies and even so it is most likely that such higher order effects will 
be overwhelmed by systematics. For the current highest reach of Run 2, E T ps 650 GeV, the 
0(a\ctw) effects studied here amount to only —2.5% for ji = E T /2. If one adopts another 
factorisation/renormalisation scale (see Fig. 18), they can be slightly larger. They increase 
in excess of —3% for /j, = Et/4. For a fixed scale, e.g., \i = Mz, they reach —4%. Notice 
however that for higher jet transverse energies the scale dependence reduces considerably. Of 
more relevance numerically is the case when the LO term is defined to be only of 0(al). As 
discussed in [23], for certain choices of /i, the cumulative effect of EW corrections through 
(9(as«EW + «ew + c^s^w) could well reach the —10 to —15% level for Et values which will 
be measurable within Run 2, with the tree-level and one-loop EW terms being of similar size. 
Thus, it is of paramount importance to establish which terms are included in Monte Carlo 
(MC) programs used to interpolate the data. 

Clearly, in the case of Tevatron, the aforementioned logarithmic effects bring very little 
enhancement, as the Et values that can be probed at such a machine are not much larger 
than M\y and Mz- Therefore, despite through 0{a\a^j) there are many more diagrams 
available for channels yielding two jets in the final state than via 0(a^) or indeed C(as«Ew) 
and C(a| w ) 6 , the fact that the Sudakov regime is not reached here implies that the size of 
the corrections to the LO terms of 0(a\) is never much larger than 0(a^) (i.e., with no 
logarithmic enhancement). Besides, at FNAL energies, jet production is always dominated 
by quark-antiquark induced channels, so that large corrections typical of processes via quark- 

6 Notice that subprocesses (5), (6), (9) and (10) are CKM suppressed in our diagonal approximation through 
O(asaiEw)- 
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quark/antiquark-antiquark scattering (see later on) do not contribute much to the hadronic 
cross section. 

5 LHC phenomenology 

More dramatic results are found at LHC energies. Here, the jet transverse energy can be 
very large so that the corrections now include the Sudhakov enhancement. Furthermore, 
there is a considerable increase in quark-quark/antiquark-antiquark luminosity which, owing 
to the large number of four-fermion Feynman diagrams through one-loop (with respect to tree- 
level), increases significantly the impact of 0(a\as N ) terms. We observe that the additional 
Feynman diagrams generally interfere constructively. For example, in our diagonal CKM 
matrix approximation, processes (9)-(10) only count one Feynman diagram at tree- level (via 
gluon exchange) whereas at one- loop they are mediated by all graphs depicted in Fig. 8. As 
some of the one-loop amplitudes are Sudakov-enhanced and some others are not and none 
of them is gauge-invariant on its own, it is not however possible to quantify a priori the 
overall impact of the increased number of graphs even channel by channel. In general, the 
larger the number of Sudakov-enhanced amplitudes the larger the corrections are. There 
is, however, a contrasting effect for the LHC due to the fact that gluon induced processes 
become dominant and these are generally subject to small weak corrections. This affects the 
overall normalisation but not the Sudakov enhancement. Altogether, despite the fact that the 
correction to the total cross section is still small, the corrections to the differential one at large 
Et are substantial, e.g., at Et = 1 TeV one-loop effects are already about —10% and increase 
further, up to —35% at 4 TeV, irrespective of the PDFs used (over most of the kinematical 
range): see Fig. 19. Also notice that the change of factorisation/renormalisation scale has a 
smaller impact at LHC than at Tevatron, as confirmed by Fig. 20. 

The shape of the corrections in the Figs. 19 and 20 can be understood in terms of the 
partonic composition of the complete 0(a$a-w) sample, distinguishing between processes ini- 
tiated by (anti) quarks only and those which have a gluon component in the initial state, see 
Fig. 21. This distinction stems from the fact that in the case of subprocesses initiated by 
(anti)quarks only, one also has LO EW effects through C(q;s«ew + «lw)- (I n the plot, the 
label LO SM identifies the sum of terms of 0(a^ + as^EW + a Ew)-) These, however, can only 
reach a 3% effect at Et = 1 TeV, with respect to the 0(a^) terms (LO QCD) and a 16% effect 
at Et = 4 TeV. On the other hand, the 0(a$ + a$aw) terms (labelled NLO weak) lead to 
corrections up to —40% in the vicinity of 4 TeV, and even at Et = 1 TeV they amount to 
—12%. Thus we observe that the one-loop corrections of order a§a: w dominate the tree-level 
interferences of order ctsctEW owing to the large logarithms. Finally, the plot also shows the 
contributions from only those subprocesses that are not initiated by any gluons (denoted by 
the label (qq)): it is clear that at very large E T (the Sudakov regime) are these channels that 
dominate much of the jet phenomenology. 

It is also of interest to understand the different behaviours of the O(ala^) effects in terms 
of each of the partonic channels involved. To this end, we present Fig. 22(a)-(c) 7 , showing 

7 Recall that individual NLO terms need not be positive definite, as they represent interferences. The total 
LO + NLO result of course is. 
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the contributions to the Et dependent cross section of subprocesses (1)— (15) (plus gg — > gg, 
qq' — > QQ', qq' — > QQ' and qq' — > at LO only) separately at the LHC. We show the 
range from E T =100 GeV, which is sufficiently far from poorly modelled threshold effects at 
Et ~ Mw/2, to _&r =800 GeV, where Sudakov effects start being active in the 0(aga^) 
corrections. The main purpose of this plot is to illustrate that, at the LHC, unlike the case 
of the Tevatron, the magnitude of the gluon luminosity inside the proton leads to gluon- 
initiated processes that can easily compete with the (anti) quark-initiated ones, even at rather 
large jet transverse energies. In fact, at LO, processes (3), (4) and gg — > gg dominate over 
the total of the (anti) quark-initiated ones for all Et values considered in the plot 8 . Tab. 1 
quantifies this, e.g., at Et = 800 GeV. The 0(a\ayj) corrections are particularly large for 
subprocesses (7)-(8), mainly by virtue of the large number of diagrams involved at loop level 
(as intimated earlier), with respect to the full LO case. Process (12) has a large positive 
correction. This can be understood from the fact that here the tree-level interference between 
strong and weak interactions is forbidden by colour and the leading contribution comes from a 
gluon exchanged in the s-channel. However, at one-loop level the cross section is substantially 
enhanced by the interference between a t-channel H^-exchange and a t-channel gluon exchange, 
with a further gluon exchanged in the t-channel in order to conserve colour. Nevertheless this 
process only makes a small contribution to the complete jet sample. Altogether, despite 
the fact that corrections to individual channels can be of a few tens of percent, the overall 
0{a\a^) effect is approximately —6%, with respect to the 0(a^asa-Ew + «ew)- Similar 
patterns for the corrections are typical also at larger jet transverse energies, with the overall 
size of the one-loop corrections increasing steadily like ~ log 2 (E T / M^) or ~ log 2 (E T /M^) 
and the quark- (anti) quark initiated processes gradually becoming dominant. 

6 Summary and conclusion 

In summary, whilst at the Tevatron 0{a\a^) one-loop virtual corrections to the single jet 
inclusive cross section may be comparable to statistical and systematic effects, at the LHC 
such terms are important contributions at large jet transverse energies. For the expected 
highest reach of the CERN machine (assuming standard luminosity), E T ~ 4 TeV, they can 
be as large as an astounding —35 to —40% (depending on whether these are related to the 
LO QCD cross section or to the total tree-level one including interference between strong and 
EW interactions and the square of the latter). Therefore, they ought to be included in any 
comparison of theory with data in these regimes. However, particular care should be paid 
to the treatment of real W- and Z-production and decay in the definition of the inclusive 
jet data sample, as this will determine whether W- and Z-bremsstrahlung effects have to be 
included in the theoretical predictions through 0{a\a^), which might counterbalance the 
negative effects due to the one- loop W- and Z-exchange estimated here [33]. (As these were 
not included in our calculation, the matter is currently under study [45].) 

8 The relative importance of (anti)quark-initiated processes originates from the combination of the valence 
quark luminosity, which is always large, and a large number of Feynman diagram for four-quark processes, as 
opposed to a gluon luminosity which decreases rapidly with increasing partonic energy combined with a small 
numbers of graphs with external gluons [21]. 
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Along the same lines, it should be recalled that NNLO EW terms ought to be investi- 
gated too, as it is well known from the Sudakov treatment that they may well be sizable in 
comparison to the NLO ones (see, e.g., Ref. [72]). 

Our results on weak corrections to the single jet inclusive cross section at hadron colliders 
are in line with the findings in several other hadronic processes and in various approximations 
(see Refs. [46]-[71] for an incomplete list limited to the SM). Altogether they should help to 
raise the awareness that LHC physics (primarily) is not always dominated by QCD effects, 
particularly in extreme kinematic regimes where new physics beyond the SM could manifest 
itself, possibly in observables that are parity violating, hence sensitive to genuine (i.e., of SM 
origin) EW corrections but not QCD ones. In our view, progress in evaluating higher order 
QCD effects should proceed hand-in-hand with that of assessing similar EW effects. 

7 Appendix 

7.1 Helicity amplitudes 

It is convenient to consider the virtual corrections in terms of helicity amplitudes. This 
approach has a number of advantages: 

1. Contributions from individual Feynman graphs can be added to the helicity amplitudes 
numerically. This allows flexibility for different analyses of the various terms entering 
the virtual corrections as well as a higher degree of control of the correctness of the 
results (as more internal tests can be enforced). 

2. The interference with the tree-level amplitudes can also be computed numerically. This 
avoids the cumbersome algebraic expressions that would be obtained if all possible in- 
terferences were computed analytically. 

3. For applications in which it is possible to polarise the incoming beams (e.g., at the 
Relativistic Heavy Ion Collider (RHIC) [22]), the contributions from different helicity 
combinations can be matched with the corresponding polarised PDFs. 

The formalism adopted for the case of subprocesses with external gluons, processes (l)-(4) 
in Sect. 2, has already been described in Refs. [21, 73] (see also [73, 74]). For the case of all other 
channels, processes (5)— (19) in Sect. 2 (including the bremsstrahlung contributions), we have 
found it convenient to adopt a different procedure. In order to describe this, let us consider, 
as an example, the amplitude 9 A^(s, t, Ai, A2, A3, A4), where Ai, . . . A4 are the helicities of the 
incoming and outgoing partons as indicated below. For process (13) of Sect. 2: 

?(pi,Ai) + gG» 2 ,A 2 ) -> QG» 3 ,A 3 ) + 0(P4,A 4 ), (21) 

9 Hereafter, s, t and u are the usual Mandelstam variables at partonic level, for which one has s + 1 + u = 
in the case of massless external particles. For brevity, we have removed here the 'hatted' notation s, t and u 
used elsewhere. 
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where q and Q are quarks of different generations (recall that, as mentioned elsewhere, to the 
accuracy to which we are working it is safe to neglect CKM mixing of flavours in the weak 
interactions) and Pi (i = 1, ...4) their four-momenta. 

The contribution to the amplitude for this process from any Feynman graph must be of 
the form 

A (s) (s,t, Ai, A 2 , A 3 , A 4 ) = ^(P2, A 2 )riu ( pi, Ai)u(p 3 , A 3 )r 2 V(P4, A 4 )5ai -a 2 5a 3 -a 4 , (22) 

where Ti, r 2 contain strings of 7-matrices (as well as propagators and an implied integral over 
loop momentum). They are in general tensorial with indices contracted between r\ and T 2 . 
Since we are interested in energy scales which are much larger than the mass of the 6-quark 
(and we are not considering t-quark production here), the quarks may be taken to be massless. 
In this case the strings of 7-matrices in r\, T 2 are vectors for self-energy or vertex correction 
diagrams and two-rank tensors for box diagrams. They may be reduced so that the most 
general forms are 

v(p 2 ,X 2 )F 1 u { p 1 ,X 1 ) = ai(s,t, Ai)Al(s, Ai)5 Ai ,-a 2 + h(s,t, Ai)-Bi(sAi)5 Ai ,-a 2 , (23) 

with 

Ai(s,Ai) = u(p2,A 2 )7-tt;i«(Pi,Ai) = (24) 

and 

£i(s,Ai) = v(p 2 , A 2 )7 • nu(p!, Ai) = iXi^/s, (25) 

where 

< = -^=(up1 + tp% + sp%) (26) 

is a unit vector in the scattering plane orthogonal to p\ and p 2 and n M is a unit vector normal 
to the scattering plane. 

The coefficients ai(s,t, Ai) and b\(s,t, Ai) have the same tensorial structure as I\ and can 
be obtained by taking traces with the appropriate projection operators, i.e., 

ai(M,A0 = ^Tr^r l7 ^i (1 ~ 2 Al75) ), 

bifrtM = ^Tr^r l7 -n (1 " 2 Al75) ). (27) 

The general form for u(p3, A3)r 2 t>(p 4 , A4) is obtained similarly 

u(p 3 , A 3 )r 2 W(p 4 , A 4 ) = a 2 {s,t, X 3 )A 2 {s, A 3 )5 A3 -a 4 + b 2 (s,t, X 1 )B 2 (sX 3 )Sx 3 -a 4 , (28) 

where 

^2(5, A 3 ) = u(p 3 , A 3 )7 • w 2 v(Pa, A 4 ) = -y/s (29) 

and 

B 2 (s, A 3 ) = u(p 3 , A 3 )7 • nv(p 4 , A 4 ) = iAsv^, (30) 
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with 



w% = —= (tp£ + up% + sp%) . (31) 
v stu 



The contribution to the complete helicity amplitude from any graph can therefore be 
specified in terms of the coefficients d\, b\, a 2 , b 2 . 

For example, consider the tree-level graph contributing to process (13) in eq. (21) due to 
the exchange of a Z-boson, 




In Feynman gauge, we set Ti = 7^ and T 2 = 7^, and multiply the entire graph by 

1 z z 
= 9qM9QM (s- M 2 z + iT z M z y 

where g%\ x 9q\ z are the couplings of the Z to the quarks q and Q, respectively. We find 

a\ = Fwi, 
bi = Fn", 
a 2 = Fw 2fl , 

b 2 = Frifj,, (32) 
so that the contribution to the helicity matrix element is 

A = F S « + ^)(- W2 , t + ?v) = 9 Z qM 9Q M1 A/r2 ~ -p rrr [s + 2t + sX^} . (33) 

The large set of Feynman graphs which contribute to any elementary process can be 
classified in terms of a small number of prototype graphs. The prototype graphs for process 
(21) are listed in the next section along with their contributions to the helicity matrix element. 

For the case of quark-antiquark annihilation in which the quarks q and Q are of the same 
flavour (or at least the same generation), there will also be contributions to the helicity matrix 
element from graphs involving the exchange of gauge bosons in t-channel. The contribution 
from such graphs can be obtained by crossing symmetry from the graphs with the gauge 
bosons exchanged in s-channel, 

A it \s,t,X 1 ,X 2 ,X 3 X A ) = -A^(t,s,X 1 ,-X 3 ,X 2 -X 4 ). (34) 

Similarly, once we have the helicity matrix element for quark-antiquark annihilation, the 
amplitudes for quark-quark or antiquark-antiquark scattering can be obtained using the usual 
crossing relations. 
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7.2 Individual topologies 



In this subsection we display the contributions to the total helicity amplitude of process (13) 
in eq. (21). The graphs are calculated in Feynman gauge in the DR scheme. Both IR and 
UV poles have been subtracted using a common subtraction scale, \x. A generic mass M is 
here attributed (for illustration purposes only) to all massive internal gauge bosons. Finally, 
all coupling constants and colour factors have been set to unity (we indicate the appropriate 
colour factor in cases where there may be a sign ambiguity). Notice that the forthcoming 
contributions have the same overall phase as the tree-level amplitude given in eq. (33). 

7.2.1 Self-energy graphs 

The contribution to the amplitude from a graph with a self-energy correction due to a massive 
gauge boson (mass M) on an external leg 



is 




The contribution to the amplitude from a graph with a gluon loop on the internal gluon line 



is 




(Here, the relevant colour factor is Ca-) 

The contribution to the amplitude from a graph with a fermion-loop self-energy correction 
to an internal gluon line 
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is 




(Here, the relevant colour factor is T R .) 
7.2.2 Vertex graphs 

The contribution to the amplitude from a correction to a vertex due to a massless gauge boson 



is 




The contribution to the amplitude from a correction to a gluon vertex due to a massless 
gauge boson with triple gauge-boson coupling 



is 




The convention for the triple gauge-boson coupling has been chosen such that the colour factor 
associated with this graph is Ca- 
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The contribution to the amplitude from a correction to a QCD vertex due to a massive 
gauge boson 



<3(A„) 



IS 



16tt 2 



where 



x 



x 




— s 



- 1 



In 



Li 2 (x) = dilog(l — x) 



2\2 



(s + M ) 



Ti- 
ll 



Li 2 1 + 



M 2 



ln(l - y) 

y 



dy. 



(40) 



(41) 



We have assumed here that all the quark masses are negligible. In the case where the outgoing 
quarks are 6-quarks and the gauge boson in the loop is a W though, the internal quark is a 
t-quark, whose mass must be taken into consideration. In such cases the contribution to the 
amplitude is not readily expressed as an analytic function and the contribution is calculated 
numerically using Veltman-Passarino reduction [75] and the FF library [76] . 



7.2.3 Box graphs 

The contribution to the amplitude from the box graph with two massless gauge bosons ex- 
changed in s-channel 



Q(A 3 ) 



\QSIQSIS1SLQJ- 



1S 



16tt 2 



x 



4vH 



+ (1+AiAg 



,S 2 , A i„ l S \ j n 



7r + 4 In 

u \t 



2 M t 



(42) 
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The contribution to the amplitude from the box graph with a massless and a massive 
gauge boson exchanged in s-channel 



Q(-v,; 



WVWA- 

is 




Note that the terms containing u in the denominator are always multiplied by a coefficient 
that vanishes when s = —t, so that there is no bogus singularity in the backward direction. 
The contributions from box graphs in which the two gauge bosons are crossed in s-channel may 
be obtained from the above box contributions by using the substitutions t «-> u and A3 — > —A3. 
In the end, all virtual corrections to all partonic processes with external quarks and/or anti- 
quarks can be obtained from these prototype graphs and expressions by multiplying by the 
appropriate couplings and colour factors and by using the appropriate crossing relations. For 
reason of space we do not report here all such results. Rather, we make available upon request 
the codes that implement them. 
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Figure 1: The set of interferences that contributes to process 1 (gg 
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Figure 2: Continuing the set of interferences that contribute to process 1 {gg 
reverse the direction of time, 2 (qq — > gg). 22 
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Figure 3: The set of interferences that contribute to 
direction of the fermion arrow, 4 (qg — > qg). 
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Figure 4: Continuing the set of interferences that contribute to process 3 (qg — > qg) and, if we 
reverse the direction of the fermion arrow, 4 (qg — > 5(7). 
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Figure 5: The set of interferences that contributes to process 5 (qq — > qq) and 6 (qq — > gg), the 
latter obtained by reversing the arrows on all fermion lines of the former. 
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Figure 6: Continuing the set of interferences that contributes to process 5 (qq — > qq) and 6 (qq — ► qq), 
the latter obtained by reversing the arrows on all fermion lines of the former. 
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Figure 7: The set of interferences that contributes to process 7 (qQ — ► (same generation)) and 
8 (qQ — ► g<5 (same generation)), the latter obtained by reversing the arrows on all fermion lines of 
the former. Here, the weak couplings to the two fermion lines will be different, this means we cannot 
implement the last two interferences as a factor of 2 as we could in Figs. 1,2. 
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Figure 8: The set of interferences that contributes to process 9 (qQ — > (different generation)) 
and 10 (gQ — > gQ (different generation)), the latter obtained by reversing the arrows on all fermion 
lines of the former. Here, the weak couplings to the two fermion lines will be different, this means 
we cannot implement the last two interferences as a factor of 2 as we could in Figs. 1, 2. 
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Figure 11: The set of interferences that contributes to process 11 (qq — > qq): continued from Figs. 9 
and 10. 
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Figure 14: The set of interferences that contributes to process 14 (qQ — > qQ (same generation)). 
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Figure 15: The set of interferences that contributes to process 15 (qQ — > qQ (different generation)). 




Figure 16: The interferences that contribute to process 16 (diagrams with a Z exchange, qg — ► qqq) 
and 18 (diagrams with a exchange, qg — ► qQQ). If we reverse the direction of the incoming 
fermion line we also obtain process 17 (qg — ► ggg) and 19 (g^ — > qQQ). 
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Figure 17: The effects of the 0(a§aw) corrections [bottom] relative to the full LO results 
(i.e., through 0(a\ + agctEw + a lw)) ft°p] f° r the case °f Tevatron (Run 2) for three choices 
of PDFs. They are plotted as function of the jet transverse energy E T . The cut \r}\ < 2.5 
has been enforced, alongside the standard jet cone requirement AR > 0.7. The factorisa- 
tion/renormalisation scale adopted was ji = jj, F = ji R = E T /2. 
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Figure 18: The effects of the O(a^a^) corrections [bottom] relative to the full LO results 
(i.e., through 0{pt\ + as^EW + a Ew)) [t°p] f° r the case of Tevatron (Run 2) for three choices 
of factorisation/renormalisation scale. They are plotted as function of the jet transverse 
energy Et- The cut \rj\ < 2.5 has been enforced, alongside the standard jet cone requirement 
AR > 0.7. The PDFs used were CTEQ6L1. 
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Figure 19: The effects of the 0{a\a^) corrections [bottom] relative to the full LO results (i.e., 
through C(«g + «s«ew + ^Iw)) ft°p] f° r ^ ne case °f LHC for three choices of PDFs. They 
are plotted as function of the jet transverse energy E T . The cut \r]\ < 2.5 has been enforced, 
alongside the standard jet cone requirement AR > 0.7. The factorisation/renormalisation 
scale adopted was /i = /ip = fiR — E T /2. 



38 



jet — production (|^| < 2.5) 
5 = (NLO-LO)/LO 




10 



-20 



-30 



1000 



2000 
E T (GeV) 



3000 



4000 



Figure 20: The effects of the O(a^aw) corrections [bottom] relative to the full LO results 
(i.e., through 0(a$ + asctEW + ct|w)) [t°p] f° r the case of LHC for three choices of factorisa- 
tion/renormalisation scale. They are plotted as function of the jet transverse energy Et- The 
cut \r]\ < 2.5 has been enforced, alongside the standard jet cone requirement AR > 0.7. The 
PDFs used were CTEQ6L1. 
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\ Q Q ^•~"~T~~ | | | | | | | | | | | | | | | | | | 

' 1000 2000 3000 4000 

E T (GeV) 

Figure 21: Top: The total single jet inclusive distribution in transverse energy through 0(a§) 
at the LHC. Bottom: The effects of the one-loop 0(a|aw) an d tree-level C(as«EW + a lw) 
corrections relative to the spectrum above. The label (qq) refers to the case of subprocesses 
with no gluons in the initial state. Rates are plotted as function of the jet transverse energy E?- 
The cut \tj\ < 2.5 has been enforced, alongside the standard jet cone requirement AR > 0.7. 
The PDFs used were CTEQ6L1 whilst the factorisation/renormalisation scale adopted was 
H = fi F = H R = Et/2. 
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Figure 22: The total single jet inclusive distribution in transverse energy through tree-level 
(LO) via 0{a\), 0(asctEw) and C(«ew) (top) and one-loop (NLO) corrections at 0{a\a^) 
at the LHC The cut \q\ < 2.5 has been enforced, alongside the standard jet cone requirement 
AR > 0.7. The PDFs used were CTEQ6L1 whilst the factorisation/renormalisation scale 
adopted was \i = fip = [Mr = E T /2. Absolute rates are given for each of the suprocesses 
(1)— (15), labelled as follows (round brakets imply summation over corresponding channels): 



(a) 


99 - 


-> qq (solid black, LO & NLO) 


(a) 


qq - 


-> g g (fine — dashed black, LO & NLO) 


(a) 


Q(Q)9 ~ 


-> q(q)g (dotted black, LO & NLO) 


(a) 


99 ~ 


-> gg (dashed black, LO only) 


(b) 


qq - 


-> qq (solid red, LO & NLO) 


(b) 


qq - 


-> qq (fine — dashed red, LO & NLO) 


(b) 


qQ - 


-> qQ (same generation) (dotted red, LO & NLO) 


(b) 


qQ 


-> gQ (same generation) (dashed red, LO & NLO) 


(b) 


qQ 


-> gQ (different generation) (dot — dashed red, LO & NLO) 


(b) 


qQ 


-> gQ (different generation) (wide — dotted red, LO & NLO) 


(c) 


qq - 


-> gg (solid blue, LO & NLO) 


(c) 


qq - 


-> QQ (same generation) (fine — dashed blue, LO & NLO) 


(c) 


qq - 


-> QQ (different generation) (dotted blue, LO & NLO) 


(c) 


qQ 


-> gQ (same generation) (dashed blue, LO & NLO) 


(c) 


qQ - 


-> gQ (different generation) (dot — dashed blue, LO & NLO) 


(c) 


q(q)q'(q') - 


-> Q(Q)Q'(Q') (wide - dotted blue, LO only) 


(c) 


qq' - 


-> QQ' (solid purple, LO only) 
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y/j = 14 TeV, E T = 800 GeV 



Subprocess 


(a) 


(b) 


(c) 


QQ — > QQ 






14.3 


(1) 


-0.0292 


-3.88 


0.643 


(2) 


-0.0189 


-3.97 


0.408 


(3)-(4) 

V / V / 


-1.089 


-1.75 


53.2 


(5)-(6) 


-1.078 


-9.36 


9.85 


(7)-(8) 


-2.78 


-23.9 


9.93 


(9)-(10) 


-0.239 


-6.68 


3.07 


(11) 


-0.169 


-5.82 


2.48 


(12) 


0.0390 


33.9 


0.0985 


(13) 


-0.0274 


-8.09 


0.290 


(14) 


-0.0861 


-2.98 


2.47 


(15) 


-0.239 


-6.69 


3.05 


qq' -> QQ' or qq' -> QQ' 






0.0616 


gg' QQ' 






0.00239 


Total 


-5.71 







Table 1: The contributions of subprocesses (1)— (15) to the total correction through 0(ag«w) with 
respect to the full LO result, i.e., through 0(a| + as«EW + a| w ), in the case of the differential 
cross section at LHC for Er = 800 GeV. Column (a) indicates the percentage contribution of the 
correction to the total; column (b) indicates the percentage correction to the individual partonic 
process; column (c) indicates the percentage contribution from that partonic process at tree-level to 
the differential cross section. Here, we have paired together the channels with identical Feynman 
diagram topology. The cut \q\ < 2.5 has been enforced, alongside the standard jet cone requirement 
AR > 0.7. The PDFs used were CTEQ6L1 whilst the factorisation/renormalisation scale adopted 
was (i = (if = fJ-R = Et/2. 
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